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WEPS 模型 在 乌 兰 布 和 沙漠 油 沙 豆 ( Cyperus 
esculentus ) 种 植 区 的 应 用 
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摘 E: 为 探究 WEPS 模 型 的 风蚀 预测 效果 ,在 内 蒙古 乌 兰 布 和 沙漠 油 莎 豆 ( Cyperus esculentus ) 种 植 区 ,对 全 采 收 


(CK)、 留 4 采 6、 留 6 采 6 油 莎 豆 纯 作 和 油 莎 豆 - 梭 梭 间 作 


留 茬 4 种 不 同 采 收 模式 地 表 进 行 风蚀 监测 ,并 利用 2020 年 


11 月 15 日 和 12 月 26 日 2 次 实测 风蚀 数据 对 WEPS 模 型 风蚀 预测 结果 进行 验证 。 结 果 表 明 :(1) 与 全 采 收 相 比 ,3 种 


留 茬 模式 均 能 有 效 降 低 风 力 对 种 植 区 地 表 的 侵蚀 作用 ,地 于 


防风 固沙 能 力 随 留 茬 数量 的 增加 而 提高 。(2) WEPS 模 


型 预测 单 宽 风蚀 量 与 实测 结果 有 明显 差异 ,模型 值 最 大 为 实测 值 的 10.16 倍 , 最 小 为 实测 值 的 0.58 倍 ;模型 预测 质 


量 有 不 确定 性 ,在 植被 覆盖 度 较 高 的 地 表 预 测 效 果 较 差 。(3) 模型 预测 的 单 冤 输 沙 量 随 实测 值 的 增加 哇 窜 函数 增 


加 ,能 够 较 准 确 地 预测 不 同 地 表 特 征 风蚀 量变 化 趋势 ;在 风蚀 定量 化 估算 上 ,WEPS 模 型 还 需 根据 实际 风蚀 环境 对 


公式 及 参数 进行 修正 。 


Xu: WEPS 模 型 ; WHYS (Cyperus esculentus); 风力 侵蚀 ; 乌 兰 布 和 沙漠 


风蚀 是 气流 (风力 ) 作 用 下 土壤 圈 或 岩石 圈 的 
破损 ,其 过 程 是 风力 作用 引起 的 地 表 物 质 脱 离 地 
表 、 搬 运 和 再 堆积 过 程 的 统一 趾 , 是 沙 质 葛 漠 化 的 主 
要 成 因 , 也 是 干旱 、 半 干旱 及 部 分 半 湿 润 地 区 的 主 
要 环境 问题 之 一 外。 准确 预测 风蚀 及 风蚀 引起 的 环 
境 变 化 ,对 于 土壤 保持 规划 指导 与 检验 各 种 土壤 
风蚀 防治 措施 减轻 风蚀 引起 的 空气 污染 .维护 风 
蚀 土 地 的 可 持续 利用 等 十 分 必要 ”。 有 目前 ,土壤 风蚀 
估算 主要 有 野外 直接 观测 所 风 洞 模拟 试验 5 、 同 位 
素 分 析 法 ”粒度 对 比 法 上 和 风蚀 模型 模拟 ”2 等 方 
法 。 通 常 认为 ,一 个 成 功 的 土壤 风蚀 模型 能 够 预测 
不 同 尺 度 和 不 同 地 表 类 型 的 风蚀 情况 ,是 估算 风蚀 
量 最 科学 dec HR SUUS 

20 世 纪 30 年 代 末 期 ,现代 流体 力学 的 创立 使 风 
蚀 的 定量 化 研究 成 为 可 能 ;随后 几 十 年 ,许多 国家 
相继 建立 起 适用 于 不 同 尺度 的 定量 模型 来 估算 风 
刨 强度 '”。 小 尺度 模型 适用 于 面积 分 布 在 100~ 
1000 的 区 域 , 具 有 代表 性 的 有 风蚀 方程 (WEQ)、 
修正 风蚀 方程 (RWEQ) 风蚀 预报 系统 (WEPS) 、 德 


收 稿 日 期 : 2022-03-04; ”修订 日 期 : 2022-05-28 


克 萨 斯 风蚀 分 析 模 型 (TEAM) 、 风 蚀 随 机 仿真 模型 
(WESS) SEU ;中 等 尺度 模型 适用 面积 最 大 可 达 
10000 km ,如 欧洲 轻 质 土壤 风蚀 模型 (WEELS) 、 风 
蚀 评 价 模型 (WEAM) 、 综 合 风蚀 模型 系统 
(IWEMS) ;大 尺度 模型 适用 于 面积 在 10000 km 以 
上 区 域 ,目前 最 重要 的 大 尺度 模型 为 粉尘 释放 模型 
(DPM)"。 其 中 ,WEPS 模 型 由 美国 农业 部 组 织 

发 ,是 目前 最 完整 和 最 先进 的 土壤 风蚀 预报 模型 ， 
它 适用 范围 广 , 不 仅 针对 农田 区 域 还 兼顾 草原 地 区 
的 风蚀 预测 所 。 想 要 大 范围 的 运用 WEPS 模型 ,还 
需要 拥有 适合 国家 级 的 WEPS 数 据 库 "5 。 美 国 农业 
部 于 2016 年 将 WEPS 软件 更 新 到 WEPS_1.5 版 本 ， 
建立 起 相对 完备 的 WEPS 数 据 库 , 使 WEPS 模 型 得 
到 广泛 运用 并 成 为 农田 风 刨 的 主要 评价 工具 。 陈 
莉 等 "在 2012 年 首次 使 用 WEPS 模 型 对 天 津 郊区 
土壤 风蚀 起 尘 量 及 迁移 量 进 行 估算 , 探 明了 郊区 起 
侍 迁 移 情 况 对 市 区 的 影响 ; 王 燕 等 "在 民 勤 荡 漠 地 
区 运用 WEPS 模型 对 风蚀 量 进行 预测 , 发 现 模型 在 
模拟 无 植被 覆盖 区 域 时 与 实测 值 比较 接近 。 刘 现 
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A5 55 FH RWEQ 模型 和 WEPS 模 型 对 北方 农 牧 交错 
带 潜在 风蚀 进行 评估 ,得 到 WEPS 模型 模拟 效果 要 
PEF RWEQ 模型。Liu 等 运用 WEPS 模型 对 中 国 
新 疆 南 部 风蚀 源 PM 和 了 PM:。 排 放量 进行 模拟 ,其 准 
确 性 较 高 。 

WH Yb (Cyperus esculentus ) 原 产地 为 非洲 干旱 
沙漠 区 ,根系 发 达 , 分 药 能 力 强 ,具有 抗 逆 性 强 、 耐 
次 薄 、 适 应 性 强 等 优异 特性 ; 除 此 之 外 ,还 具有 防风 
固沙 及 地 力 培育 等 生态 修复 功能 。 通 过 科学 种 植 
油 莎 豆 及 合理 构建 油 莎 豆 防风 固沙 采 收 模式 ,能 
效 遏 制 北方 风沙 区 土壤 自然 侵蚀 和 表层 沙土 流动 ， 
为 解决 我 国 北方 沙漠 化 土地 风蚀 沙化 严重 ,水 资源 
短缺 . 盐 渍 化 加 剧 等 问题 提供 可 能 2。 有 效 的 风蚀 
估算 是 构建 油 莎 豆 科 学 种 植 模式 的 前 提 , 但 传统 凡 \ 
蚀 监 测 手 段 , 因 工 作 量 大 、 监 测 时 段 长 和 监测 范围 
有 限 等 特点 ,加 大 了 风蚀 监测 工作 的 难度 。 研 究 
WEPS 模 型 在 内 蒙古 油水 豆 种 植 区 的 运用 效果 不 仅 
为 该 区 域 提 供 一 种 高 效 的 风蚀 监测 方法 ,也 将 为 我 
国运 用 风蚀 模型 在 小 尺度 区 域 进行 风蚀 预测 提供 
理论 基础 。 


1 研究 区 概况 


PRR MPA SH AIA Ee sR TE RE HE 
(图 1) ,地 处 马兰 布 和 沙漠 东 缘 ,土地 面积 4167 
km’, 沙漠 面积 2847 km 。 海 拔 在 1030-2046 m, 除 
北部 狠 山 山脉 外 ,地势 整 体 由 东南 向 西北 降低 ,地 
SAVE .平原 和 山地 为 主 ;研究 区 属 温 带 大 陆 性 
季风 气候 ,年 平均 气温 7.6 % ,年 平均 降水 量 145 
mm ,年 平均 蒸发 量 2398 mm ,水 热 时 空 分 布 不 均 ,年 
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图 1 研究 区 位 置 示意 图 


Fig. 1 Sketch map of study area location 
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际 变化 较 大 。 土 壤 质 地 较 轻 ,土壤 类 型 包括 风沙 
ERE KEE cf T WEE UB DR 
BY VAUD 3-293. X P HEL BC F BL. TRAN RUBRIC TRIER 
(Salix chaenomeloides ) 胡杨 (Populus euphratica ) Jy 
FR (Ulmus pumila ) FE (Haloxylon ammodendron ) , 
等 ;草本 与 灌木 植物 有 白 刺 (WNitraria tangutorum ) 、 
沙 打 旺 (hstragalus adsurgens) , Z T E& ll (Tamarix 
ramosissima ) Fr 2& 6X5 JL (Caragana korshinskii) 、 细 
AX 4 BF (Corethrodendron scoparium ) FU W W (Hip- 
pophae rhamnoides ) ^5 o 


2 材料 与 方法 


2.1 试验 设计 

在 乌 兰 布 和 沙漠 油 莎 豆 种 植 示 范 区 内 选择 3 个 
油 莎 豆 纯 作 地 块 及 1 个 油 莎 豆 和 梭 梭 间 作 地 块 开 展 
人 研究 ,4 个 地 块 相 邻 且 独 立 , 大 小 均 为 80 mx140 m, 
地 块 内 油 莎 豆 种 植 方 式 为 带 状 种 植 ,每 人 巷 2 行 油 莎 
豆 , 纯 作 种 植 方向 为 东 一 西 走向 ,间作 呈 南 一 北 走 
向 。 将 1 个 纯 作 地 块 内 的 油 莎 豆 进 行 全 部 采 收 以 作 
对 照 ( 全 采 收 ,CK) ,其 余地 块 均 在 中 部 建立 了 不 同 
采 收 方式 的 试验 小 区 ,小 区 内 的 采 收 类 型 分 别 为 : 
EI 4 EK 628 ARO) E16 227622 C EH 626) AL 
2E lb EA BO 1 TRIR BER MCCIRWE RI TE)3 fj 
采 收 模式 ,大 小 分 别 为 34 mx40 m、39 mx40 m,20 
mX40 m。 并 于 2020 年 11 月 15 日 和 12 月 26 日 ,对 4 
种 采 收 模式 (图 2) 地 表 分 别 进行 了 2 次 风蚀 观测 ,各 
模式 内 呈 阶 梯 式 分 别 均匀 布设 3 个 集 沙 仪 ( 集 沙 口 
大 小 为 1 cmx3 em FAI D 4 16 AH EE ) RH HO-BO 
多 剖面 自 计 式 风 速 仪 (监测 距 地 面 高 度 0.5 m, 1 m 
1.4m 处 风速 ); 每 个 模式 内 选取 4 个 1mx1 m 样 方 ， 
测定 植被 直立 株 高 . 叶 面 积 指数 带宽 , 带 间 距 \ 植 
株 密度 ,并 在 样 方 内 采用 五 点 交叉 法 进行 土壤 取 
样 ,分 别 采用 环 刀 法 和 烘 干 法 测定 0~15 cem EE 
壤 容 重 及 土壤 含水 率 。 在 下 午 13:30 一 14:00 在 人 研究 
区 上 空 用 无 人 机 拍摄 光学 正 射 影像 ,拍摄 高 度 10 m, 
焦距 为 9 mm, 最 大 光圈 2.97 , 图像 空间 分 辨 率 为 7 
mm; 利用 Agisoft PhotoScan 软件 进行 影像 拼接 处 
HE , TE ENVIS.3 中 采用 颜色 混合 分 析 法 (CMA) ZR 
取 植 被 盖 度 。WEPS 模型 中 结 皮 面 积 比 表示 为 除 植 
被 外 的 地 表 可 见 覆 盖 物 使 地 表 不 起 沙 面积 比 , 人 研究 
中 将 覆盖 在 各 小 区 地 表 的 可 见 滴灌 带 面 积 测算 为 
结 皮 面 积 ,通过 在 ENVI5.3 中 对 无 人 机 影像 采用 支 
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图 2 采 收 模式 布局 及 试验 布设 示意 图 


Fig.2 Schematic diagram of harvesting mode and instrument layout 


持 向 量 机 法 进行 监督 分 类 计算 得 到 。2 次 风蚀 事件 
中 ,不 同 采 收 模式 风速 情况 及 地 表 特 征 参 数 分 别 见 
表 1 和 表 2。 
2.2 WEPS 模型 

风蚀 预报 系统 (Wind Erosion Prediction System, 
WEPS ) 为 模块 化 结构 ,由 1 个 用 户 界面 .1 个 管理 程 
序 7 个 子 模型 和 4 个 数据 库 组 成 ,是 基于 过 程 的 每 
日 时 间 步 长 模型 ,不 仅 可 以 模拟 基本 的 风蚀 过 程 ， 
还 可 以 模拟 改变 土壤 对 风蚀 敏感 性 的 田间 管理 和 
风化 过 程 。 目 前 ,WEPS_1.5 版 本 无 中 国土 壤 、 作 物 
和 耕作 管理 方式 等 基础 数据 文件 ,不 能 在 软件 中 直 


接 选 择 参数 模拟 风蚀 过 程 。 本 文采 用 WEPS 模 型 中 
的 风蚀 子 模型 对 日 风蚀 量 进行 计算 ,首先 确定 阔 值 
摩 阻 风速 (模型 中 最 小 值 设 为 0.35 m-s”) , 当 风 速 超 
过 阔 值 时 ,计算 代表 该 地 表 的 一 系列 单个 网 格 单元 
的 土壤 风蚀 量 。 

风蚀 量 预 测 公式 

Q,, = 0.4U.(U. - 0.8U.,) (1) 

AP: 0, 为 单 宽 输 沙 量 (kg.m');U, 为 摩 阻 风速 
(m-s); Us, 为 临界 摩 阻 风速 ms )。 

摩 阻 风速 U, 的 计算 : 

(1) 计算 气象 站 摩 阻 风速 U. 


表 1 不 同 采 收 模式 2 次 风速 监测 


Tab.1 Twice wind speed monitoring under different harvesting modes 


日 期 /年 -月 -日 ” 采 收 模式 ”0.5 m 人 处 日 平均 风速 /(m:s?) 1.0m 处 日 平均 风速 /ms …) 1.4 mAb HSEJS G/(m s). 起 风 时 长 h UY 
2020-11-15 ”全 采 收 6.88 7.83 7.87 22 西 偏 北 
2020-12-26 6.09 6.85 7.23 23 
2020-11-15 留 4 采 6 6.78 7.29 8.00 22 西 偏 北 
2020-12-26 5.33 6.30 6.85 23 
2020-11-15 ” 留 6 采 6 6.37 7.03 7.12 22 西 偏 北 
2020-12-26 3.93 4.78 5.89 23 
2020-11-15 间作 留 茬 6.50 7.04 7.36 22 西 偏 北 
2020-12-26 3.81 4.95 571 23 
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表 2 不 同 采 收 模式 地 表 特 征 


Tab.2 Soil surface characteristics of different harvesting modes 


采 收 模式 平均 密度 /( 株 "m”) 带宽 /m 带 间 距 /m 土壤 容重 /(g'em”) EREKE% «50 km 颗粒 /% 
全 采 收 0 0 1.64+0.02a 1.78+0.47c 25.72+0.52a 
留 4 采 6 52.00+3.00a 4.20+0.23a 5.50+0.16a 1.57+0.04a 2.34+0.28b 16.17+0.34¢ 
E636 48.00+2.00b 5.50£0.35b 5.50+0.21a 1.51+0.01b 2.24+0.43b 24.71+0.41a 
间作 留 茬 47.00+3.00b 0.70+0.06c 2.50+0.11b 1.60+0.02a 2.9420.53a 19.59+0.46b 


_ 0.4U 
In 


of 
式 中 : U 为 高 度 为 Z 时 的 实测 风速 ,为 10 m fe BE JA 
3E (m « s) ; Zy 为 气象 站 空气 动力 学 粗糙 度 ,在 
WEPS 中 取 25 mm; 
(2) 计算 观测 地 点 摩 阻 风速 
D 计算 无 植被 覆盖 时 , 摩 阻 风速 


Z 0.067 


式 中 : Z 为 当地 空气 动力 学 粗糙 度 ,在 本 文 分别 对 
应 4 种 采 收 模式 近 地 表 的 空气 动力 学 粗糙 度 。 

Q) 计算 有 植被 覆盖 时 , 摩 阻 风速 

当地 表 有 直立 植被 时 , 需 考 虑 直立 植被 对 摩 阻 
风速 的 前 弱 作 用 ,通过 判断 有 效 植 被 拖 忠 系数 的 大 
小 选择 合理 的 公式 进行 计算 。 

(a) 有 效 植被 拖 忠 系数 

BR = BR,,(0.2 -0. 15e ^) +BR,, — (4) 


AP: BR, 为 有 效 生物 量 拖 忠 力 系数 ; BR 为 直立 
植被 叶 面积 指数 (mm?); BR. 为 茎 面积 指数 (mx… 
m?) ,表示 茎 面积 与 水 平地 表面 积 的 比值 , 油 莎 豆 是 
草本 植物 ,在 此 忽略 不 计 。 

(b) 植被 冠 层 空 气动 力学 粗糙 度 

根据 有 效 植被 拖 忠 系数 是 否 大 于 0.1, 判 断 直立 
植被 冠 层 空气 动力 学 粗 燃 度 计算 公式 如 下 。 


BZ 
Z = , (BR,,>0.1 
0 A 1.254 In(BR.,) a Hey) 
. B BRA 7 BR, 
(5a) 


BZ 4.60517 0.001 
Z,= (Sb) 


(0.001 « BR,, <0.1) 
0 , (BR,, <0.001) 


0.11 a 
Z Her B 
BZ 0 + BZ w Ra ， 


注 : 同 列 不 同 小 写字 母 表示 不 同 采 收 模式 间 差 异 显 著 (P<0.05)。 下 同 。 


WP: Zo 为 直立 植被 冠 层 空气 动力 学 粗糙 度 (m ) ; 
BZ 为 直立 植被 高 度 (m)。 
(e) 植被 冠 层 上 部 摩 阻 风速 
Z 0.067 
U =U (6) 
: (2) 


(d) 植被 冠 层 下 部 摩 阻 风 速 


0:250, e% 


, BR, »2.56BZ 


U,- -Bu -BR (7) 
2 * msc BR. <2.56BZ 


临界 摩 阻 风速 0 的 计算 : 
(1) 光滑 平坦 地 表 摩 阻 风速 
UB, -1.7-135e " 
a 1 
WT (8) 


式 中 : SF, 为 地 表 履 盖 土 块 / 结 皮 或 石 块 的 面积 百 
分 比 。 
(2) 地 表 有 倒 放 植物 引起 的 临界 起 动 摩 阻 风 速 


SC, = BF,,(1-SF,,) 
UC,,, = 0.02 + SC,, 
式 中 : SC, 为 倒 放 植被 引起 的 临界 起 动 摩 阻 风速 的 
增 量 系数 ;BF 为 倒 放 植被 覆盖 率 ; UC. 为 倒 放 植 
被 引起 的 临界 起 动 摩 阻 风速 增 量 (ms )。 
(3) 含水 率 引 起 的 临界 起 动 摩 阻 风速 增加 量 


H, 


How H 2 
ov. -asi -1 -o4 owe J (10) 


(9) 


HR. 


RP: HV UHR KA (kg kg '); HR 为 1.5 
MPa 地 表 含 水 率 (kg*kg')。 
(4) 总 临界 起 动 摩 阻 风速 
U., = UB, + UC., + UW, (11) 
WEPS 中 临界 起 动 摩 阻 风速 的 最 小 值 设 为 0.35 
m*s"'。 因 此 ,可 根据 计算 所 得 摩 阻 风速 与 总 临界 起 
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动 摩 阻 风 速 的 关系 ,采用 式 (1) 计 算 单 宽 输 沙 量 。 
3 结果 与 分 析 


3.1 实测 风蚀 结果 

3.1.1 输 沙 量 随 高 度 分 布 规律 ”从 2 次 风蚀 事件 的 
风沙 流 结构 特征 (图 3) 来 看 , 随 着 油 莎 豆 留 荐 数量 
的 增加 ,风沙 流 结构 特征 由 曲折 状 逐 渐变 为 “1” 字 
形 。 第 1 次 风蚀 事件 中 ,在 0~0.21 m 输 沙 高 度 内 ,全 
采 收 模式 地 表 不 同 高 度 的 输 沙 量 均 高 于 其 他 3 种 留 
茬 模式 ;但 在 0.21~0.39 m 内 ,4 种 模式 地 表 不 同 高 度 
的 输 沙 量 差异 较 小 , 且 随 高 度 增 加 ,各 层 输 沙 量 趋 
近 于 0( 图 3a)。 与 全 采 收 模式 相 比 ,3 种 留 茬 模式 均 
能 有 效 降低 0.21 mm 以 下 的 地 表 输 沙 量 ,表现 出 较 好 
的 固沙 效果 。 第 2 次 风蚀 事件 的 风沙 流 结构 特征 与 
第 1 次 相似 ,但 3 种 留 茬 模式 的 输 沙 高 度 由 0.21 m 
降低 至 0.15m 以 下 (图 3b)。 结 合 2 次 风蚀 事件 来 


42 
39 Læ (a) 2020-11-15 一 一 AI 
36 $ —e- 间作 留 茬 
33 到 - 留 4 采 6 
30 —— 留 6 采 6 
27 

E 24 1 
21 

TE 18 + 
15 Lé» 
12 AA 
9 ÜL à 
6L4 Ba CSS C 
3L u œi BR 
0 L L 1 L L 

0 1 2 3 4 5 


HW RI (g-cm?-d^) 


看 , 随 着 留 茬 数量 的 增加 ,不 同 高 度 处 的 输 沙 量 呈 
现 :全 采 收 > 留 4 采 6> 留 6 采 6> 间 作 留 茬 的 趋势 (图 
3)。 结 果 表 明 ,间作 留 茬 模式 的 固沙 能 力 最 强 , 增 
加 留 荐 数量 可 以 提高 固沙 能 力 。 

从 不 同 高 度 与 输 沙 量 的 最 优 拟 合 关系 来 看 ( 表 
3) ,在 2 次 风蚀 事件 中 ,全 采 收 . 留 4 采 6 和 留 6 采 6 
模式 地 表 的 输 沙 量 随 高 度 分 布 均 呈 指数 关系 , 且 全 
采 收 模式 的 拟 合 度 最 高 ,其 次 为 留 4 采 6 模 式 , 最 后 
为 留 6 采 6 模式 ;间作 留 茬 模式 呈 对 数 规律 ,第 2 次 
风蚀 事件 的 拟 合 度 比 第 1 次 低 。 表 明 在 留 荐 区域， 
输 沙 量 随 高 度 呈 指数 分 布 ,但 随 着 留 荐 数量 的 增 
加 , 拟 合 度 逐渐 降低 ,并 向 对 数 函 数 方向 变化 。 
3.1.2 不 同 采 收 模式 单 宽 输 沙 量 单 宽 输 沙 量 为 某 
一 时 段 内 通过 单位 宽度 的 总 输 沙 量 , 可 以 表征 风蚀 
强度 ,该 值 越 大 ,地 表 所 受 风蚀 作用 越 强 。 通 过 计 
算 1d 内 4 种 不 同 采 收 模式 0~0.39 m 高 度 范围 内 的 


39 |a (b) 2020-12-26 —e- 全 采 收 
36 Le —e- WERTE 
33 Le —a— 留 4 采 6 
30 1 一 一 留 6 采 6 
27 上 
24 » 
21 上 
HE 18 a 
15- . 
12 上 e 
N Ky 
: L Taji inc" 
0 L 1 1 L L | 
0 1 2l 3 4 5 


Eg cmd’) 


图 3 不 同 采 收 模式 输 沙 量 随 高 度 分 布 特征 


Fig. 3 Sediment discharge distribution rules with height under different harvesting modes 


表 3 不 同 采 收 模式 输 沙 量 与 高 度 的 拟 合 关系 


Tab.3 Fitting relationship between sediment transport and height in different harvesting modes 


日 期 /年 -月 -日 采 收 模式 拟 合 方程 R P 
2020-11-15 全 采 收 Q =6.3343e °?” 0.970 4.4E-3 
留 4 采 6 Q=9.342e 7 0.947 6.2E-3 
留 6 采 6 Q=0.5458e *^* 0.923 4.6E-3 
间作 留 茬 Q 70.0567 - 0.017 In(Z) 0.920 7.8E-5 
2020-12-26 全 采 收 Q=12.919e 2 0.978 1.4E-3 
留 4 采 6 Q-1.1503e "^ 0.948 22E-2 
Fa 6X6 ( 2 0.876e * 7 0.921 5.3E-3 
间作 留 茬 Q 70.1376 - 0.039 In(Z) 0.907 2.7E-5 
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Qs= 2.0, (12) 


RE: 0, 为 单 宽 输 沙 量 (g.cm!…d'); Q, 为 第 m 个 
梯度 输 沙 量 (g.cm'…d!)。 

经 式 (12) 计 算得 到 4 种 不 同 采 收 模式 下 地 表单 
宽 输 沙 量 ( 表 4)。 在 2 次 风蚀 事件 中 ,不 同 采 收 模式 
间 的 地 表单 宽 输 沙 量 有 显著 差异 ,全 采 收 模式 单 宽 
输 沙 量 最 大 ,间作 留 荐 模式 单 宽 输 沙 量 最 小 ;在 第 1 
次 风蚀 事件 中 ,与 全 采 收 模式 相 比 , 留 4 采 6、 留 6 采 
6 和 间作 留 茬 模式 分 别 减 少 了 74.61% , 84.5096 和 
98.76% 的 单 宽 输 沙 量 ; 第 2 次 风蚀 事件 中 , 留 4 采 6、 
留 6 采 6 和 间作 留 茬 分 别 减 少 了 69.64% , 84.8696 f 


表 4 不 同 采 收 模式 下 地 表单 宽 输 沙 量 
Tab.4 Unit-width sediment discharge of different 


harvesting modes /(g* em '*d') 


日 期 /年 -月 -日 ”全 采 收 留 4 采 6 留 6 采 6 PMR 


2020-11-15  10.06x0.02a 2.55+0.05b 1.56+0.05c 0.13+0.01d 


2020-12-26 10.43+0.10a 3.17+0.03b 1.58+0.04c 0.36+0.02d 
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96.57% ,与 第 1 次 风蚀 事件 相 比 ,间作 留 茬 和 留 4 采 
6 的 固沙 能 力 轻 微 减 弱 。 
3.2 WEPS 模型 预测 结果 
3.2.1 模型 输入 预测 参数 2 次 风蚀 监测 均 在 油 莎 
豆 枯萎 期 ,无 灌溉 ,天 气 晴朗 ,气温 较 高 ,人 研究 在 使 
用 WEPS 模 型 时 ,忽略 土壤 水 分 的 影响 。 根 据 3 个 
风速 梯度 数据 ,用 对 数 廊 线 法 计算 不 同 采 收 模式 地 
表 的 空气 动力 学 粗糙 度 。10m 高 度 处 的 风速 通过 
全 采 收 模式 地 表 观 测 到 的 风速 值 拟 合 风速 廓 线 得 
到 , 因 4 种 模式 处 于 同一 研究 区 ,10 m 高 度 处 的 风速 
几乎 不 受 地 表 植 被 的 影响 , 故 取 值 一 致 。 模 型 中 输 
入 不 同 采 收 模式 地 表 的 风蚀 预测 参数 见 表 5。 
3.2.2 模型 风蚀 预测 ”通过 将 各 采 收 模式 地 表 所 测 
的 空气 动力 学 粗糙 度 、 叶 面积 指数 .植被 高 度 等 指 
标 带 入 式 (1)~ 式 (11) ,计算 得 到 2 次 风蚀 事件 中 ,4 
种 采 收 模式 的 摩 阻 风 速 ` 有 效 植被 拖 搜 系数 和 植被 
冠 层 空气 动力 学 粗糙 度 等 预测 参数 ( 表 6)。 

模型 中 的 摩 阻 风速 是 指 上 层 风 速 经 过 植被 及 
其 障碍 物 削 弱 后 作用 于 地 表 的 风速 ,能 直观 反映 植 


表 5 模型 中 输入 参数 


Tab. 5 Parameters about soil surface used in WEPS 


ar a) aT E ea euR 
2020-11-15 全 采 收 11.40 0.004 0 0 0 0 0 
2020-12-26 9.40 0.002 0 0 0 0 0 
2020-11-15 留 4 采 6 11.40 0.032 0.844 0.196 18.003 1.561 0 
2020-12-26 9.40 0.013 0.624 0.154 15.323 2.177 0 
2020-11-15 留 6 采 6 11.40 0.062 0.672 0.153 28.832 1.728 0 
2020-12-26 9.40 0.047 0.432 0.132 24.126 1.328 0 
2020-11-15 间作 留 茬 11.40 0.097 3.055 0.251 48.429 1.673 0 
1.974 0.218 45.362 1.824 0.12 


2020-12-26 9.40 0.062 


表 6 模型 参数 输出 结果 
Tab.6 Parameters calculated with WEPS 


日 期 /年 -月 -日 采 收 模式 摩 阻 风速 /m's" ) 有效 植 被 拖 忠 系数 植被 冠 层 空 气动 力学 粗糙 度 /m ”临界 摩 阻 风速 /(ms"') 


2020-11-15 全 采 收 0.673 0 0.370 
2020-12-26 0.530 0 0.370 
2020-11-15 留 4 采 6 0.739 0.177 0.020 0.673 
2020-12-26 0.615 0.124 0.018 0.642 
2020-11-15 E636 0.728 0.134 0.018 0.809 
2020-12-26 0.601 0.084 0.013 0.682 
2020-11-15 间作 留 茬 0.734 0.611 0.021 0.884 


2020-12-26 0.689 0.395 0.020 0.730 
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被 及 其 障碍 物 在 垂直 方向 上 对 风速 的 削弱 能 力 ,其 
值 越 小 表示 植被 及 障碍 物 对 风速 的 削弱 能 力 越 
强 。 模 型 计算 各 采 收 模式 下 摩 阻 风速 大 小 表现 为 : 

间作 留 茬 > 留 4 采 6> 留 6 采 6> 全 采 收 ,其 中 间作 留 茬 
模式 摩 阻 风速 最 大 ,主要 是 因为 油 莎 豆 - 梭 梭 间 作 


。 综 合 分 析 模 型 各 项 风蚀 预测 参数 得 到 ,3 种 留 
-o ppl auf 间作 留 荐 模式 
防风 效果 最 优 , 留 4 采 6 防风 效果 最 差 。 

3.3 实测 单 宽 输 沙 量 与 模型 预测 值 
从 4 种 采 收 模式 地 表单 宽 输 沙 量 的 预测 情况 


留 荐 模式 的 植被 结构 较 均 匀 , 琉 透 性 好 , 穿 过 冠 层 
内 部 的 气流 较 多 所 导致 。 在 3 种 植被 覆盖 模式 中 ， 
间作 留 茬 有 效 植被 拖 搜 系数 和 植被 冠 层 空气 动力 
学 粗糙 度 最 大 , 留 6 采 6 最 小 。 

WEPS 模 型 预测 的 临界 摩 阻 风速 是 地 表 沙 粒 开 


(图 4) 可 以 看 出 ,不 同 采 收 模式 的 单 宽 输 沙 量 实测 
值 与 模型 值 差 异 较 大 ;其 中 间作 留 荐 模式 地 表 差 异 
最 大 ,模型 预测 结果 最 大 为 实测 值 的 10.16 倍 ,最 小 

为 实测 值 的 0.58 倍 。 就 模型 预测 效果 而 言 ,第 1 次 
风蚀 事件 中 ,模型 对 全 采 收 模式 地 表 预 测 最 优 ; 在 


始 移动 的 风速 浆 值 。 在 2 次 风蚀 预测 中 ,临界 摩 阻 “第 2 次 风蚀 事件 中 ,模型 低估 了 全 采 收 模式 地 表 的 
风速 大 小 均 为 :间作 留 茬 > 留 6 采 6> 留 4 采 6> 全 采 单 宽 输 沙 量 ,对 留 4 采 6 模式 地 表 预 测 最 优 ;但 2 次 
16 12 

一 一 模型 值 

14 Li = 实测 值 10 
ae 2 -- 模型 值 /实测 值 
T 8 
2 K 
ins) 8 6 地 
e 1 / ,8 
a 4 E TS E 

; - E |, 

全 采 收 Woko 留 6 采 6 T 全 采 收 HOR Moko BEME 
2020-11-15 2020-12-26 


bas] 同 采 收 模式 


4 单 宽 输 沙 量 模型 值 与 实测 值 对 比 


Fig.4 Comparison ofmodel and measured results of unit-width sediment discharge 


风蚀 事件 中 ,模型 对 间作 留 茬 模式 地 表 预 测 均 为 最 

差 ,说 明 WEPS 模型 不 适用 于 植被 覆盖 度 较 高 的 地 
表 风 蚀 预 测 。 虽 然 ,预测 效果 在 2 次 风蚀 事件 中 有 
轻微 差异 ,但 模型 预测 地 表单 宽 输 沙 量 结果 与 实测 
基本 同步 变化 , 均 呈 现 全 采 收 > 留 4 采 6> 留 6 采 6> 间 
作 留 茬 的 趋势 。 分 析 实 测 值 与 模型 值 的 拟 合 关 系 
(图 $) ,采用 决定 系数 ( 尼 ) 和 均 方 根 误差 (RMSE ) 对 
模型 精度 进行 综合 评价 。 从 图 5 中 可 以 看 出 ,模型 
值 随 实测 值 的 增加 呈 震 函数 增加 ,2 次 风蚀 事件 尼 
在 0.91 以 上 ,模型 对 第 1 次 风蚀 事件 预测 的 均 方 根 
误差 高 于 第 2 次 风蚀 事件 。 表 明 WEPS 模 型 预测 结 
果 与 实测 结果 具有 较 高 的 相关 性 ,模型 对 第 2 次 风 
刨 事件 的 预测 精度 高 于 第 1 次 。 


^ 2020-11-15 
e 2020-12-26 
一 一 乘 短 (2020-11-15) 


? a c 乘客 (2020-12-26) y=4.06x058 
g e 
E! 10 RMSE=4.51 
m 8 
gi 
ak 6 
4 y=1.8312x949 
R?-0.98 
2 RMSE-2.25 
0 1 1 | 1 r 
2 4 6 8 10 12 
Scil f/(g-cm" d") 
图 5 单 宽 输 沙 量 模型 值 与 实测 值 函 数 关系 
Fig. 5 Functional relationship between model and measured 


value of unit-width sediment discharge 


202212.00151v1 


chinaXiv 


4 讨论 


根据 马兰 布 和 沙漠 油 莎 豆 种植 区 4 种 不 同 采 收 
模式 地 表 实 际 发 生 的 风蚀 过 程 , 对 WEPS 风蚀 模型 
预测 效果 进行 验证 。 单 宽 输 沙 量 的 模型 值 与 实测 
值 存在 一 致 性 ,这 种 一 致 性 是 非 线 性 的 ,与 学 者 们 
在 德国 和 美国 运用 风蚀 预测 系统 对 农田 风蚀 量 模 
拟 研 究 结论 基本 一 致 ”“”。 模 型 不 仅 能 够 有 效 预 测 
土壤 质地 较为 复杂 的 农田 地 表 风 蚀 量变 化 趋势 ,而 
且 对 以 沙 粒 为 主 的 油水 豆 种 植 区 地 表 也 同样 具有 
较 好 的 模拟 效果 。 同 时 ,模型 预测 的 单 宽 风蚀 量 随 
地 表 植被 盖 度 的 增加 而 降低 ,与 学 者 们 利用 WEPS 
模型 模拟 不 同情 景 地 表 下 风蚀 量 的 研究 相 吻合 ” ， 
模型 能 够 较为 真实 地 反映 油 莎 豆 种 植 区 风力 对 不 
同 植被 盖 度 下 地 表 物 质 的 侵蚀 机 制 ,也 能 够 较 准 确 
地 预测 马兰 布 和 沙漠 油 莎 豆 种 植 区 不 同 采 收 方式 
地 表 的 固沙 能 

从 模型 对 地 表 风 蚀 量 预测 精度 来 看 ,模型 预测 
质量 具有 不 确定 性 。 首 先 ,在 同 1 次 风蚀 事件 中 LBS 
型 预测 的 单 宽 输 沙 量 与 实测 结果 具有 明显 差异 ;其 
次 ,在 2 次 风蚀 事件 中 ,模型 预测 不 同 采 收 模式 地 表 
风蚀 量 的 精度 差异 较 大 。 一 方面 可 能 是 对 风蚀 事 
件 的 观测 次 数 不 足 , 以 及 在 同 1 次 风蚀 事件 中 对 不 
同 采 收 模式 地 表 的 风蚀 观测 重复 较 少 ,从 而 引起 实 
测 数 据 不 稳定 导致 ; 男 一 方面 ,由 于 阵风 的 风向 、 风 
速 和 风 频 的 不 稳定 性 影响 实际 获取 的 风蚀 量 与 真 
实情 况 不 符 汪 ,所 产生 的 监测 误差 导致 ; 且 WEPS 建 
模 实 施 每 小 时 步 长 ,并 未 考虑 阵风 的 侵蚀 效应 ,可 
能 造成 模型 预测 结果 偏 低 ™; 此 外 ,相对 于 模型 建 
立 在 植被 空间 结构 均匀 的 背景 而 言 , 人 研究 区 地 表 哇 
密集 条 状 的 植被 结构 能 有 效 防止 “ 狭 管 效应 ”的 发 
生 ”, 从 而 减少 气流 穿 过 植被 冠 层 筷 际 时 ,流速 增 
大 ,加 速 对 地 表 侵 蚀 的 现象 出 现 , 这 也 是 导致 模型 
预测 单 宽 输 沙 量 偏 高 的 重要 原因 之 一 。 

综合 来 看 ,WEPS 模 型 预测 的 单 宽 风蚀 量 对 风 
速 `. 地表 和 植被 特征 响应 较为 敏感 。 模 型 在 小 尺度 
预测 不 同 区 域 的 风蚀 量变 化 趋势 上 是 可 行 的 ,在 指 
导 农 田 风 蚀 管理 中 具有 重要 作用 。 而 在 风蚀 定量 
化 研究 上 ,WEPS 模 型 预测 精度 不 佳 , 可 通过 缩小 风 
蚀 预测 时 段 .增加 观测 频次 和 选择 持续 大 风 天 气 等 
方法 来 提高 模型 预测 精度 ,但 这 些 方法 都 是 基于 改 
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变 外 部 环境 来 更 好 地 响应 模型 ,对 模型 的 运用 和 发 
展 具 有 一 定局 限 性 ;未 来 要 真正 提高 WEPS 的 普 适 
性 ,还 需要 加 强 对 模型 在 多 气候 、 地 表 植 被 覆盖 多 
类 型 ,预测 多 时 序 以 及 其 他 与 风蚀 相关 因素 中 的 研 
究 ,并 建立 数据 库 , 从 而 根据 实际 风蚀 环境 对 模型 
参数 及 公式 进行 修正 。 


5 结论 


(1) 与 全 采 收 相 比 ,3 种 留 荐 模式 均 能 有 效 降低 
风力 对 种 植 区 地 表 的 侵蚀 作用 ,其 中 ,间作 留 茬 防 
风 固 沙 能 力 最 强 , 留 6 采 6 最 弱 , 地 表 防 风 固沙 能 力 
随 留 茬 数量 的 增加 而 提高 。 

(2) 全 采 收 、 留 4 采 6 和 留 6 采 6 输 沙 量 随 高 度 
增加 呈 指 数 降低 的 变化 趋势 ,间作 留 茬 输 沙 量 随 高 
度 增加 呈 对 数 降低 的 分 布 规律 ; 随 留 茬 数 量 的 增加 
指数 拟 合 度 逐 渐 降 低 , 输 沙 量 随 高 度 的 分 布 规律 向 
对 数 方向 演变 。 

(3) WEPS 预测 与 实测 结果 具有 一 定 差 异 ,模型 
值 最 大 为 实测 值 的 10.16 倍 ,最 小 为 实测 值 的 0.58 
倍 ; 模 型 预测 质量 具有 不 确定 性 ,在 植被 覆盖 度 较 
高 的 地 表 预 测 效 果 较 差 ,但 存在 合理 的 一 致 性 , 模 
型 预测 的 单 宽 输 沙 量 随 实测 值 的 增加 明 短 函数 二 
加 。 模 型 能 够 较为 准确 地 预测 不 同 地 表 特 征 的 风 
蚀 量 变化 趋势 ,对 于 风蚀 定量 化 合算 而 言 ,WEPS 模 
型 还 需 建立 数据 库 ,根据 实际 风蚀 环境 对 公式 及 参 
数 进行 修正 。 
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Application of the Wind Erosion Prediction System in the Ulan Buh Desert 


Cyperus esculentus planting area 


LIAO Guiyun, WU Xiuqin, TAN Jin, LI Dan, FENG Mengxin 
(School of Soil and Water Conservation, Beijing Forestry University, Beijing 100083, China) 


Abstract: The purpose of this study was to explore the wind erosion evaluation accuracy of the Wind Erosion 
Prediction System (WEPS). Wind erosion monitoring was carried out on the surface of four different harvest 
modes in the Ulan Buh Desert Cyperus esculentus planting area in Inner Mongolia, China as follows: full harvest, 
retain four harvest six, retain six harvest six, and retain intercrop (Cyperus esculentus and Haloxylon 
ammodendron). Wind erosion monitoring was carried out on the model surface, and the wind erosion prediction 
results of the WEPS model were verified using the wind erosion data, which was measured twice. The research 
results showed that compared with full harvesting, the three stubble modes effectively reduced the effect of wind 
erosion on the surface of the planting area. The intercropping stubble had the strongest wind- proof and sand- 
fixing ability among them, and the six remaining six harvested was the weakest. The surface wind-proof and sand- 
fixing ability was the weakest. With the increased number of remaining stubble, the structural characteristics of 
wind-sand flow gradually changed from a tortuous shape to the shape of *1". This can provide a basis for the 
promotion of planting Cyperus esculentus in northern sandy areas to balance ecological and economic benefits. 
With the increased height, the sediment load of full harvest, with four remaining six harvested and six remaining 
six harvested, showed an exponential decreasing trend. The intercropping stubble sediment load showed a 
logarithmic decrease with the increase of height. The fitting degree decreased gradually, and the distribution law 
of sediment transport with height changed to a logarithmic direction. There was a certain difference between the 
WEPS prediction and the measured results. The maximum model value was 10.16 times the measured value, and 
the minimum was 0.58 times the measured value. The model prediction quality was uncertain, and the prediction 
effect was poor on the surface with high vegetation coverage. However, there was reasonable consistency as the 
unit-width sediment load predicted by the model increased as a power function with the increase of the measured 
value. Overall, the model is feasible in small-scale areas to predict the variation trend of wind erosion with different 
surface characteristics. However, in the quantitative prediction of wind erosion, the prediction value is less stable. 
It is necessary to strengthen the research on the model in multi-climate, multi-type surface vegetation coverage, 
multi-sequence prediction, and other factors related to wind erosion to establish a database that adjusts the model 
parameters and formulas according to the actual wind erosion environment to improve the universality of WEPS 
in the future. 


Keywords: Wind Erosion Prediction System; Cyperus esculentus; wind erosion; Ulan Buh desert 


